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ارزیابی ریسک رو گذری سد تبار کآباد بر پایه تحلیل فراوانی یک متغیره و دو متغیره سیلاب* 
آتوسا رضازاده عنبرانی ٩‏ محمود فغفور مغربی ۲ 

چکیده از میان دلایل مهم شکست سد حاکی روگذری عمده‌ترین دلیل است. در تحلیل شکست سد به علت روگذری تحلیل عدم قطعیت و 
محاسبه ریسک ضروری است. در این پژوه شکاربرد تحلیل ریسک شکست سد حاکی بر اثر سرریز شدن آب از بدنه سد بر پایه تحلیل فراوانی 
یک متغیره و دو متفیره سیلاب با استفاده از توزیع لوگ رمال برای سد اک تبارکآباد واقع در خحراسان رضوی بررسی شده است. ننایج تحلیل 
فراوانی دو متغیره سیلاب در شش هیدروگراف جریان ورودی و با در نظ رگرفتن 00۷ (ترکیب دبی و حجم سیلاب) برای دورة بازگشت‌های 
توأم ۰۶ ۰ و ۲۰۰ سال ارائه شده است. ریسک روگذری بر اساس تحلیل فراوانی تک متغیره سیلاب و برای هم هیدر وگراف‌های حاصل از 
تحلیل فراوانی دو متغیره و شش تراز اولی آب در محزن با درنظرگرفتن دبی اوج سیل, تراز اولیه آب و ضریب دبی سرریز به‌عنوان متغیرهای 
غیرقطعی ارزیابی شده است. برای انجام تحلیل عدم قطعیت از روش شبیه سازی مون تکارلو و تکنیک نمونهگیری لائین هایپ رکیوب استفاده شده 
است. نهاییأ نتایج به‌دست‌آمده با استفاده از هر دو روش تحلیل فراوانی سیلاب یک متغیره و دومتغیره در دورة بازگشت‌های مختلف با هم مفایسه 
گردید. نتایج ذشان داد که ریسک روگذری برآورد شده در همه دور با زگشت‌ها بر اساس تحلیل فراوانی سیلاب دو متغیره بیشتر از روش تک 
متغیره بوده که عطرات بیشتری به همراه دارد. همچنین هیدروگراف ۷1-۵0 که دارای بیشترین حجم سیلاب بوده, در همهٌ ارتفاع‌ها از ریسک 
بیشتری دسبت به بقیه هید روگراف‌ها برخوردار است. در برآورد ریسک روگذری بر اساس تحلیل فراوانی یک متغیره و دو متغیره», بالارفتن تراز 
آب در مخزن, نسبت به استفاده از سیلاب‌های با دورة بازگشت طولانی‌تر منجر به روند افزایش ریسک محسوس‌تری حواهد شد. به‌طور یکه در 
ارتفاع 11 متری تراز سط حآب» مقدار ریسک روگذری برای سد بسیار بالا است. 

واژه‌های کلیدی عدم قطعیت. ریسک روگذری, تحلیل فراوانی یک متغیره و دو متغیره سیلاب. شبیه‌سازی مونت‌کارلو نمونه‌گیری لاتین 
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(۱) دانش أموختة کارشناسی ارشد آب و سازه‌های هیدرولیکی, گروه مهندسی عمران, دانشکده مهندسی. دانشگاه فردوسی مشهد 
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۸۴ 


مقدمه 

در سدها یا سازه‌های کنترل سیلاب برای جلوگیری از 
هرگونه فروپاشی و يا شکست ناگهانی سازه باید توجه 
ویژه‌ای به دبی سیلاب شود. به‌عنوان مثال با طراحی مناسب 
سرریز سد و ظرفیت کنترل سیلاب یک مخزن می‌توان ایمنی 
سد را اطمینان بخشید و از هر گونه مشکلات نامطلوب مانند 
روگذری جلوگیری کرد. از این‌رو برآورد دقیق دبی طراحی 
سیلاب و هیدروگراف جریان برای طراحی سازه‌های 
هیدرولیکی مورد نیاز است. دبی طراحی سیلاب به‌عنوان 
حداکثر جریان سیل که یک سازه می‌تواند با اطمینان عبور 
دهد تعریف می‌شود [1]. در ارزیابی ریسک سیلاب از آنالیز 
تک متغیره استفاده می‌شود که تنها متغیر دبی اوج سیلاب 
مدنظر قرار می‌گیرد اما این روش یک تکنیک خیلی دقیق 
نیست و موجب دستیایی به اطلاعات محدود در زمینه‌ی 
ریسک سیلاب می‌گردد. درحالی‌که در برنامه‌ریزی و 
طراحی‌های هیدرولوژیکی مرتبط با مدیریت سیلاب تنها 
آگاهی نسبت به متغیر دبی اوج سیلاب کافی نیست و نیاز به 
اطلاعات و آگاهی در زمینه سایر متغیرهای مهم یعنی تداوم 
و حجم سیلاب نیز است. برای دسترسی به اطلاعات در 
زمینه متغیرهای سیلاب با استفاده از تحلیل‌های آماری چند 
متغیره. نیاز به توابع توزیع تجمعی و توابع چگالی احتمال 
توأم متغیرهای سیلاب است [۲]. 

بسیاری از تلاش‌ها برای تحلیل فراوانی دو متغیره و 
چند متغیره سیلاب با در نظر گرفتن ارتباط بین متغیرهای 
مختلف سیلاب از جمله دبی اوج. حجم سیلاب و 
مدت‌زمان سیل انجام شده است. عصقصصتت [3] وعع00ظ و 
«عوصصصعم؟۳ [4] مطالعات جامعی درباره‌ی استفاده از تحلیل 
فراوانی سیلاب تک متغیره ارائه دادند و به این نتیجه رسیدند 
که تحلیل فراوانی تک متغیره نمی‌تواند ارزیابی دقیقی از 
پارامترهای سیل ارائه دهد. از این‌رو باید تحلیل فراوانی دو 
متغیره يا چند متغیره سیلاب برای مشخص کردن بهتر 
هیدروگراف جریان با در نظر گرفتن سایر عوامل مهم مانند 
حجم رواناب و یا مدت‌زمان سیل به کار گرفته شود. 6061 
وهمکاران اوق اواتش تسیادت اه شیر بو که ارت اند 


نشريه مهندسی عمران فردوسی 


ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک... 


از حجم سیلاب. دبی اوج و تداوم جریان با استفاده از توزیع 
نرمال تحلیل کردند. ۷۷۵ [6] از توزیع نرمال دو متغیره برای 
تحلیل توزیع احتمال همبسته توأّم دبی پیک و حجم و 
توزیع احتمال همبسته توأم حجم و مدت‌زمان سیل استفاده 
کرد. ۷۷۵ [7] با استفاده از توزیع لوگ نرمال دو متغیره 
روشی برای توزیع احتمال همبسته توأم دبی پیک و حجم 
و توزیع احتمال همبسته توأم حجم و مدت‌زمان ارائه کرد. 
همچنین توزیع‌های توأم توزیع‌های شرطی و دور 
بازگشت‌های مرتبط با این متغیرهای تصادفی از توزیع‌های 
حاشیه‌ای به دست آمد. 1۷/0616 6 و همکاران [8] یک 
توزیع احتمالی دو متغیره را با استفاده از تابع مفصل و توزیع 
مقدار حدی تعمیم‌یافته دو متغیره برای دو پارامتر دبی سیل 
و حجم ارائه کردند. مفصل‌ها توابعی هستند که توزیع‌های 
چند متغیره را به توزیع‌های حاشیه‌ای یک متغیره پیوند می- 
دهند. همچنین, رفتار مخزن با استفاده از تعداد زیادی از 
هیدروگراف‌های سیل برای بررسی ایمنی هیدرولوژیکی 
مخازن و کفایت سرریز سد در نظر گرفته شده بود. 
مقصصه ۷ و تصنصت [9] اطمینان‌پذیری روگذری 
یک سد را با استفاده از تحلیل فراوانی سیلاب دو متغیره 
ارزیابی کردند. تحلیل فراوانی دو متغیره سیلاب با در نظر 
گرفتن ترکیب دبی پیک سیلاب و حجم رواناب و برای 
دور بازگشت ۱۰۰ سال و با استفاده از توزیع گاما دو متغیره 
انجام شد. از طریق مطالعٌ آن‌هاء حداکثر ارتفاع مخزن و 
ریسک روگذری با انجام احتمالاتی روند یابی مخزن بر 
اساس روش شبیه‌سازی مونت‌کارلو تعیین شد. منصوری و 
همکاران [۱۰] به‌منظور توسعه کاربرد آنالیز ریسک در 
تحلیل ایمنی سدها؛ به محاسبهٌ خطر روگذری سد خاکی 
ونک با در نظر گرفتن عدم قطعیت پارامترهای موثر در این 
پدیده با استفاده از روش‌های شبیه‌سازی مونت‌کارلو و لاتين 
هایپرکیوب پرداختند. نتایج نشان داد که در نظر گرفتن عدم 
قطعیت این پارامترها در احتمال شکست سد موثر بوده 
است. همچنین دقت در برآورد بزرگی بی پیک سیلاب نقش 
بسزایی در کاهش ریسک روگذری سد دارد. درحالی که 
تأثیر تراز اولیه آب بر میزان خطرپذیری سد قابل اغماض 
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آتوسا رضازاده عنبرالی- محمود فغفور مغربی 


است. بر اساس ننایج آنالیز ریسک با تلفیق اثر باد و سیل 
میزان ریسک سد به‌طور محسوسی افزايش می‌یابد. گودرزی 
و همکازان [11] با استفاده از تعلیل:فراوانی, سلات لک 
متغیره و توزیع احتمالاتی گامبل لاجستیک دو متغیره برای 
دور بازگشت ۱۰۰ سال و با استفاده از شبیه‌سازی مونت- 
کارلو و لاتين هایپرکیوب. احتمال روگذری سد درودزن در 
جنوب ایران را بررسی کردند. نتایج نشان داد که ریسک 
روگذری با افزایش سطح اولیه آب در هر دو روش عدم 
قطعیت افزایش می‌یابد. ذهبیون و شرافتی [۱۲] ریسک 
روگذری سد را با درنظرگرفتن عدم قطعیت بارش و سیلاب 
تحلیل کردند. آنها با استفاده از مدل‌های هن۳) 13۳0 
(0606۲2۱00) معالد۲ و 1۳01 سیلاب با در نظر گرفتن 
عدم قطعیت‌های بارش و پارامترهای مدل بارش - رواناب 
تولید نمودند. همچنین با استفاده از شبیه‌سازی مونت‌کارلو 
تحلیل ریسک روگذری سد جامیشان در استان کرمانشاه انجام 
شده که اثرات ابعاد سررین ارتفاع سد و تراز نرمال بر ریسک 
بررسی شده است. نتایج نشان داد که در مقایسه دو عامل 
هیدرولیکی مثر بر سرریز شدن سد. حجم کنترل سیلاب بر 
تک کشت | باتش کت خر هراق 

سالاری و هم‌کاران [۲] از توابع مفصل برای 
مدل‌بندی ساختار همبستگی و همچنین برآورد توابع توزیع 
تجمعی شرطی و نیز دوره‌های بازگشست توأم ترکیب‌های 
مختلف متغیرهای سیلاب در ایستگاه هیدرومتری اهواز بر 
روی رود خانه کارون استفاده کردند. آنها همچنین بیان 
کرد ند که طرح این سناریو های متفاوت می‌تواند برای 
ارزیابی ریسک مرتبط با مسائل هیدرولوژیکی از قبیل 
طراحی سرریز و کنترل سیلاب مفید باشد. نتایج 
ذشان‌دهندهة کارایی و انعطاف‌پذیری کاربرد توابع مفصل در 
تحلیل فراوانی دو متغیره سیلاب بود. 

همان‌طور که گفته شد بیشتر مباحث مطرح شده در 
زمینه‌ی ریسک روگذری سدهاء | ستفاده از تحلیل فراوانی 
تک متغیره مسیلاب بر اس‌اس دبی اوج برای توا ید 
هیدروگراف ورودی به‌عنوان پارامتر عدم قطعءیت بوده 
است. در اين مطالعه» ریسک روگذری و تحلیل عدم 


قطعیت بر اساس مفهوم فراوانی تک متغیره و دو متغیره 
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سیلاب برای سد تبارک‌آباد واقع در خراسان رضوی ارائه 
شده است. عوامل ا صلی غیرقطعی در تحلیل فراوانی تک 
متغیره سبلاب عبارت‌اند از: دبی اوج سیل (600) عمق اولیه 
آب در مخزن (۲۷0) و ضریب تخلیه سرریز(00) و عوامل 
غیرقطعی در تحلیل فراوانی دو متغیره سیلاب. عمق اولیه 


آب در مخزن و ضریب تخلیه سرریز می‌باشند. 


توزیع دو متغیره لوگ نرمال 
در توزیع لوگ نرمال دومتغیره. تابع چگالی احتمال ۱ 
حمناعصنط تمصع بانلزطه۳:۵0) لوگ نرمال دو متغیره 
به‌صورت زیر تعریف می‌شود: 


(3-) ««پس 


-ل(وید,0] 
۵ 27۳520۷/1 


(1+ > > 1-) 
)۱( 
که در آن: 


0 < 


با 6 
) 9 ( 1 


شرا 1۳ 
ك )+ 0 


0۷ 


2 


ربا - (م16/ بل - 1۳60 
۱ ۷2 ( ۱ 0 آّ 


2 


با و ,6۷ به ترتیب میانگین و انحراف‌معیار 
متغیرهای تصادفی در حالت لگاریتمی (11,2, 26 م[< (۷) 


که از فرمول‌های (۲) و (۲) به دست می‌آیند [13]: 


وه 
۳ پگ +1 ]| < به 
۹ 1 
۳ 9 
۳ سد - (پرلا) 1 دا 


که در آن‌ها زیلز و زره میانگین و انحراف معیار متغیر 
تصادفی :۸6 برای 12 بوده وا برایر با صریب همیسدکین 
دو متغیر تصادفی ۰۱ و ۷ است که به‌صورت زیر محاسبه 
می‌شود: 


(2 


لته 


۵2 2 


نشریه مهندسی عمران فردوسی 


و 


تابم توزیع تجمعی (6 مات :0۲ 
حمناعصنظ جمننطاتنامز) توزیع دو متغیره لوگ نرمال 
(۳06,,6 به‌صورت آنالتیک قابل حل نیست و باید 
به‌صورت انتگرال عددی محاسبه شود. دورة بازگشت برای 
دو متغیر ۱ و 262 که به ترتیب با و :1 نشان داده 


می‌شوند. می‌توان به‌صورت زیر نوشت: 


(۵) (ب > ما - (و۳0) مت < بط 
10( (ارد > ]۲ < () سس و" 


دورة بازگشت توأم به دو حالت "یا" و "و" و دور 
بازگشت شرطی با توجه به هدف طراحی تقسیم می‌شوند. 
محاسبات مربوط به دورة بازگشتی که به "یا" معروف 
است که در آن یکی از متغیرهای تصادفی ,26 یا 262 یا 
هردو از حدود آستانه :* پا 52 بیشتر شوند. که به صورت 
زیر نوشته می‌شود: 


۹2 


1 


( ارو 1-۳ 3 ود 
دوره با تین که به "و" معروف تحت تین از 
متغیرهای تصادفی :26 و 262 از حدود آستانه :و 62بیشتر 


شوند که به شکل زیر بیان می‌گردد [14]: 


1 


۳ سل و1 
1 
(0۵) ود ور (002) ور (ب0 1 و1 


(۸) 


شرایط انواع دیگر دورهٌ بازگ شت دو متغیره تو سط 
۵۵ [14] و 21720071 و ۱0616 16 [8] ارائه شده 
اش تا در این پژوهش. نوع "یا" دوره بازگشت توأم در 
تحلیل فراوانی دو متغیره سیلاب استفاده شده است. 


روند یابی مخزن 
هدف اصلی از تحلیل پدیده شکست ناشی از سرریز شدن 
آب از بدنهٌ سد خاکی؛ تخمین ارتفاع آب در مخزن با توجه 
به دبی ورودی و مقایسه با ارتفاع تاج است. در این راستاء 
سیلاب درون مخزن روند یابی می شود که ساده‌ترین فرم 


آن به‌صورت زیر نوشته می‌شود: 


نشریه مهندسی عمران فردوسی 


ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 


1۷ 


7 2100-000 (۵ 


که در آن )1 دبی ورودی به مخزن» 060 دبی 
خروجی از مخزن. ۷ حجم آب در مخزن و ] زمان است. 
می‌توان با حل معادله )٩(‏ با استفاده از روش رانگ-کوتا 
(هاانک-عط18) مرتبه چهارم حداکثر ارتفاع آب در مخزن 
را محاسبه کرد. می‌دانیم: 
(0۰ ۳7 ۶ (۸)۲ < 0۷ 


کفقر ان (۸ سطح پلان مخزن در رقوم 1 است. 
با جایگزینی معادله (۱۰) در معادلهةٌ تست کر (معادله )٩‏ 
نتیجه زیر حاصل می‌شود: 


10-0۵ _ اه 


3 ۸0۳0 (۱۱) 


با در نظر گرفتن رابطه سرریز به‌جای (06 حاصل 


می‌شود: 
10-2 _ 1۳ 
۳( مه ٩‏ 


که در آن بآ به‌عنوان طول تاج و 11 تراز سطح آب در 
دریاچه سد نسبت به تراز تاج و 00 ضریب تخلیه سرریز 
است؛ بنابراین: 


0 


11 
۱۳( (1) مه 


حال با | ستفاده از ضرایب رانگ-کوتا که به صورت 
زیر هستند و معلوم بودن :11 که ارتفاع آب در مخزن در 
گام زمانی ذام است. می‌توان :1 مخزن را در گام زمانی 
1+1 به دست آورد: 


1 
رو (یا + 216 + 2162 + و6 2+ را ع وبرلا 


که در آن: 
(رر ۵ -< 1 
(ر6ا + را رگ + ر؟< ۵۱ < و16 
روک + را رگ + نگ < ۵۱ < وکا 
(یع1 + راب۵ + 6 <۵۲ 2 16 
(۱۵( 


سال سی و دوم شمارة دو ۱۳۹۸ 


آتوسا رضازاده عنبرانی- محمود فغفور مغربی 


با معلوم شسدن با مقدار ,۵ از رابطة سبرریز 
محا سبه می شود. به همین ترتیب می‌توان 1 ها و دبی‌های 
خروجی را برای گام‌های زمانی بعدی به دست آورد [15]. 
برای افزايش دقت نتایج و همچنین کاهش اثر غیر قطعی 
بودن ] ۸ می‌توان فا صله زمانی را کاهش داد. از این‌رو در 
این مقاله. گام زمانی ۳۰ دقیقه‌ای انتخاب شده است. 


شکست یک سیستم زمانی به وقوع می‌پیوندد که سیستم 
قادر به انجام عملکرد قابل انتظار نبا شد. شکست می‌تواند 
باری (1) تعریف شود که از مقاومت و یا ظرفیت سیستم 
(8) فراتر رود. شناسایی بار وارده و مقاومت. اساس آنالیز 
ریسک بوده و به نوع ساختار هیدرولیکی و فیزیک مسئله 
بستگی دارد. احتمال شکست به‌صورت ( < )۲ 
ریسک نیز می‌تواند به شرح زیر نشان داده شود: 


0 < ۲51 < ۳0۶ > 0( (۱2۱ 


که در آن 2 تابع عملکرد است که تابعی از اجزا 
توزیع آن می‌تواند | شکال متفاوتی دا شته با شد که در زیر 
امن الست: [ 117 


7 1- (2) < 2:() ها < 2زا - 8 < 2 


مدل‌سازی ریسک برای روگذری سیل 

روگذری به حالتی گفته می‌شود که بخش خروجی سد قادر 
درگاش تاها رم ان از ره 
تاج سد بالاتر بیاید. در تجزیه‌وتحلیل روگذری. حداکثر 
ارتفاع آب در مخزن (روم3) و ارتفاع سازه‌ای سد(ع]) به 
ترتیب می‌توانند بار و مقاومت سیستم درنظرگرفته شوند. 
اگر (11) ارتفاع بالاآمده آب مخزن به‌وسیلةٌ نیروی سیل 

(۳)و ع1] ارتفاع تاج سد و وا ارتفاع سطح اولیة آب مخزن 
قبل از وقوع ۴ باشد آنگاه روگذری زمانی رخ می‌دهد 
که! < وا +11۶ شود (شکل ) [18]. لذا ریسک 


سال سی و دوم. شمارة دی ۱۳۹۸ 


۸۷ 


روگذری ۳به‌صورت زير خواهد بود: 


رم ما۲ < ۲۱۳ 2 [,۲1 < و11 +۲۱۳۲ <۲۶ 


شبکل ۱ ثمای شماتیک ر گر سل 


از این‌رو تابع عملکرد به‌صورت زیر نوشته می‌شود: 
3 
)۱۹( ‌ ) ۱ < 7 


حول 


که در آن م7تابع عملکرد سیل و مرآ آبالاترین سطح 
آب در طول یک رویداد سیل بر اساس روند یابی مخزن 


محاسبه شده است. در نهایت. ریسک روگذری به‌صورت 


زير محاسبه می‌شود: 
۳( (908 - 1 - (ک) ۵ 1 < ون 
که در آن 8 نشان‌دهنده شاخص اطمینان‌پذیری بوده و 
به‌عنوان نسبت میانگین تابع عملکرد ببلایه انحزاف معیار 
استانداردیر» تعریف می‌شود. همچنین (0)8 احتمال نرمال 
تجمعی مربوط به 8 است (شکل ۲) [18]. 


احتمال روگذری 
0-( مس )19 2 جریان ورودی 


جریان خروجی 


شاخص اطمینان پذیری (]۵ 


شکل ۲ برآورد ریسک روگذری سد [18] 


هدف از آ نالیز عدم قطء‌یت تعیین ویژگی های آ ماری 


خروجی‌های سي ستم و یا برر سی اثر عوامل تغییرپذیر در 
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عکس‌العمل سیستم است. انتخاب روش مناسب آنالیز عدم 
قطعیت به ماهیت م سئله, اطلاعات در د سترس, پیچیدگی 
مسئله و نوع نتایج مورد نظر بستگی دارد. متداول‌ترین 
روش های مورد استفاده برای آ نالیز عدم قطعءع یت در 
مدل‌های هیدرولیکی و هیدرولوژیکی روش شبیه‌سازی 


مونت‌کارلو و روش نمونه‌گیری لاتين هایپرکیوب است. 


شبیه‌سازی مونت کارلو 
(صمتاهام‌هنگ ماه عاصما0۷) 
یکی از روش‌های ارزیابی عدم قطعیت و آنالیز ریسک که 
کاربرد زیادی در علوم مختلف دارد روش شبیه‌سازی 
مونت‌کارلو است. فرآیند تولید ورودی‌های تصادفی متعدد 
برای انجام شبیه سازی‌های متعدد برای تعیین خروجی‌های 
تصادفی متعدد. شبیه‌سازی مونت‌کارلو (۷105) نامیده 
می شود. با افزايش قابل توجه در قابلیت‌های کامپیوتری و 
توسعه تکنیک‌های کاهش واریانس در سال‌های اخیره 
استفاده از این روش در زمینه‌های مختلف علمی افزایش 
یافته است. فرآیند مونت‌کارلو یک شبیه‌سازی عددی برای 
برگرداندن متغیرهای آماری بر اساس یک توزیع آماری 
ی اس دیا وی ای تا سای مش ها 
غیرقطعی موثر در عملکرد سامانه با توجه به تغییرات 
طبیعی این متغیرهاه یک توزیع احتمالاتی برای هریک در 
نظر گرفته می‌شسود. برای استفاده از هر توزیع احتمالاتی 
باید پارامترهای آن مشخص شوند. بعد از تعیین پارامترهای 
توزیع‌های مورد نظر» برای هر متغیر غیرقطعیء در هر تکرار 
از این روش یک سری اعداد تصادفی برای هریک از 
متغیرها در سامانه تولید می شود و سامانه با استفاده از هر 
ترکیب از اعداد تولید شده شبیه‌سازی می‌شود. روند اجرای 
5 زمانی به پایان می‌رسد که همگرایی لازم برای 
خروجی مدل شبیه‌سازی حاصل شده باشد. میزان همگرایی 
در این روش به تعداد تکرارها و همچنین تعداد متغیرهایی 
که به‌طور تصادفی تولید می‌شوند بستگی دارد. به‌طوری که 
هر چه تعداد متغیرها پیشتر باشد. مدل با تکرار بیشتر به 


تعواب همگر | می‌شود. همسنن با اف انشن تعداد: تکرازها 
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ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 


شرایط همگرایی بهبود می‌یابد. ولی تعداد شبیه سازی‌های 
نجام گرفته در فرآیند تحلیل افزایش می‌یابد و بدین ترتیب 
مدت ز مان اجرای فرایند 1۷05 طولانی‌تر می‌شود. از 
م ه اتقاف ی ها تاماهتا 
به‌منظور دستیابی به همگرایی. یکی از عوامل مهم در 


نمونه گیری لاتین‌هایپر کیوب 

0 0 
برای افزايش دقت خروجی‌های شبیه‌سازی مونت‌کارلوه 
روش های کاهش واریانس مختلفی وجود دارد که بدون 
افزايش حجم نمونه اين کار را انجام مید هد [1815]16 
تک او وف ماه کاهت ی وارتاشی است اه کارا 
پارامترهای آ ماری خروجی را افزایش می‌د هد. در روش 
قالط هگا تبونگیری از یک تنم # میرم ی ان 
محدوده‌ها با احتمال وقوع 1/0 صورت می‌گیرد که البته 
ها کیک همه کی دق رف رای یک 
توزیع احتمال مطلوب 1.15 یک مقدار صادفی از هریک 
از محدوده متغیر های بدون تکرار را استخراج می‌کند. 
سرعت همگرایی 1.۲15 از روش شبیه‌سازی مونت‌کارلو 
بیشتر است [17]. الگوریتم نمونه‌گیری متغیرهای تصادفی 
مستقل (0 با استفاده از تکنیک 1.215را می‌توان به‌صورت 
زیر خلاصه کرد: 
مرحله ۱: تقسیم‌بندی محدوده متغیرهای ورودی به مقدار 
1 
مرحله ۲: تولید متغیرهای یکنواخت (0,1) ۷ - ۷1 
مره ۳ میحاسبة زط از معادله زیر 


۳۱( هر. 1 6( ) ج اد ها 


مرحله : تعیین متغیرهای تصادفی د از 0۳ معکوس به 


۳» 2 ۳ ]۳[ (۳۲۲ 
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محدوده مورد مطالعه (سد تبارک آباد) 
سد مخزنی خاکی - س‌نگ‌ریزه‌ای تبارکآباد در ۲۵ 
کیلومتری شمال شرقی شهر قوچان در ۱/۵ کیلومتری 
پایین‌د ست رو ستای تبارک‌آباد که یکی از سر شاخه‌های 
فرعی رودخانه اترک است. احداث گردیده است. ساختگاه 
سد تبارک در مختصات جغرافیایی "58042730 طول 
شرقی و 3771030 عرض ش‌مالیء در هجاورت 
رو ستای تبارک آباد واقع گردیده ااست. تبارک آباد تأمین و 


تنظیم آب مورد نیاز را برای ۱1۵۰ هکتار از اراضی پایاب 


۸۹ 


و همچنین تأمین قسمتی از نیاز آب شرب درازمدت 
قوچان به میزان ۱۴ میلیون متر مکعب در سال در صورتی که 
تأمین آب شرب درازمدت شهر قوچان از منابع آب‌های 
زیرزمینی و آهکی امکان‌پذیر نباشد. فراهم می‌کند. در 
جدول (۱) خلاصه‌ای از مشخصات اصلی سد تبارک‌آباد 
نشان داده شده است» همچنین شکل (۱) بیانگر منحنی‌های 
اک باه امه 


جدول ۱ مشخصات سد تبارک آباد 


مشخصات سد توضیحات 
نوع سد خاکی سنگریزه‌ای با هسته رسی 
مساحت حوضه ۰ کیلومتر مربع 
نام رودخانه تبارک آباد 
تراز نرمال از سطح دریا ۵ متر 
تراز تاج ۵ سمتر از سطح دریا 
عرض تاج 1 
طول تاج ۸ متر 


مشخصات سد توضیحات 
حجم آب قابل تنظیم ۵ میلیون مترمکعب 
حجم مخزن در تراز نرمال | 1۰ میلیون متر مکعب 
حجم مرده ۳ میلیون متر مکعب 
نوع سرریز آزاد با تبدیل نوع بادبزنی 
ظرفیت سرریز ۷ مترمکعب بر ثأانیه 
طول تاج سرریز ۹ 
ارتفاع سد ۶ متر 


حجم مخزن (میلیون عتر عکعب) 


۲ ۲ 1 1 ۲ 
500 600 700 200 900 1000 


1525 


1920 


1515 


1910 


1505 


ارتفاع مخزن (متر) 
5 
8 


1495 


15 


1475 0 ۲ ۲ ۲ 
0 100 200 3200 4200 


سطح مخزن (هکتار) 


شکل ۳ تغییرات سطح و حجم مخزن سد تبارک‌آباد در برابر تراز ارتفاعی 
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از بین عوامل شناسایی شده اثرگذار در رفتار سیستم‌های 
آبگیری و طغیان سد. پارامترهایی چون دبی ورودی. تراز 
اولیة آب مخزن. منحنی تراز حجم مخزن سد. ظرفیت 
سرریز و... می‌توانند مشمول عدم قطعیت‌های گوناگون 
باشند. پارامترهای غیرقطعی در نظرگرفته شده در این 
مطالعه عبارت‌اند از: 
سیلاب حداکثر با دورة بازگشت‌های مختلف (م00): 
مهم‌ترین دلایل برای انتخاب این پارامتر به‌عنوان متغیر عدم 
قطعیت خطا در ثبت داده‌هاه فقدان داده و وجود دبی‌های 
جانبی ورودی به مخزن است. این متغیر تنها در تحلیل 
فراوانی سیلاب تک متغیره. غیرقطعی در نظرگرفته می‌شود. 
همچنین با استفاده از تحلیل فراوانی سیل. مقادیر ميانگین و 
انحراف معیار استا ندارد دبی اوج برای سیل با دوره 
باز گشت‌های مختلف با استفاده از نرم‌افزار ۳15 11720 
در حدول (۲) ارائه شده است. 

برای سیل با دور بازگشت‌های ۵۰ ۱۰۰ و ۲۰۰ 
سال, حداقل. حداکثر و متوسط هیدروگراف ورودی به 
مخزن سد تبارک‌آباد به ترتیب در شکل (4) نشان داده شده 


اننتنت: 


الف) (۷۲6۵۷) 1*50 


۱ 1 0 
حداقل هیدروگراف ورودی ما٩‏ 
و تون ی 13525 ۱ 500 
متوسط هیدروگراف ورودی سس 1 


حداقل هیدروگراف ورودی 
حداکثر هیدروگراف ورودی < : < 


متوسط هیدروگراف ورودی 


ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک... 


جدول ۲ میانگین و انحراف معیار 0 در دوره باز گشت‌های 


مختلف 
(9/) 0۵ 
دور بازگشت (سال) ِ_ 
میانگین انحراف معیار 
0۰ ۳/۶۶( ۸۸ 
۱.۰ 2۰/۳۸ ۱۱۵ 
۳۰۰ ۰/۹۶ ۱-۳۵۹( 


تراز اولیةٌ مسطح آب (۲۷0) : متوسط تراز سطح آب 3 
مخزن سد بر اساس سطح آب مشاهده شده و ثبت‌شده در 
مدت ۱۳ سال (۱۳۹۲-۱۳۸۰) محاسبه شده است. میانگین 
و انحراف معیار | ستاندارد عمق آب به ترتیب ۵۰/۱ متر و 
۸ متر بودند. علاوه بر اين؛ پنج عمق دیگر به‌عنوان عمق 
اولیه و به‌منظور در نظر گرفتن اثر تغییرات عمق اولیه بر 
احتمال روگذری فرض شده‌اند (جدول ۳). 


جدول ۳ مقادیر حداقل عمق اولیه آب 


٩۰ | ۶ | ۲ "۰-۰ ۱ ۸۱ 


عمق میانگین (0) 


۳ ضریب تخلیه سرریز (6/0): بر اساس گزارش‌های فنی 
سد تبارک آباد ميانگین و انحراف معیار | ستاندارد ضریب 
تخلیه به ترتیب ۲/۰۸ و ۰/۰۹۹ در نظر گرفته شده است. 


800 

ب) (۲6۵۲) 1100 ۱ 200 

حداقل هپدروگراف ورودی 312 13255 15353 5552 ۳۳ 
حداکثر هیدروگراف ورودی - ۰ سم ۱ ۱ ۱ 

متوسط هیدروگراف ورودی سس ۱ ۱ 0 

1 : 

300 

200 

س ۱ 100 


شکل ۶ حداقل. حداکثر و متوسط هیدروگراف ورودی برای دورهٌ بازگشت‌های الف) ۰ ب) ۰ وج) ۰ سال 


نشریه مهندسی عمران فردوسی 


سال سی و دوم شمارة دو ۱۳۹۸ 


(ورصب 


آتوسا رضازاده عنبرالی- محمود فغفور مغربی 


مشسخصات پارامترهای ورودی مانند میانگین () 
انحراف‌معیار استاندارد (6) و تابع توزیع احتمال (۱۲ظ) 
که منا سب برای داده ‏ صادفی غیرقطعی بودند. در جدول 


( ارات له اسر 


جدول ۶ خصوصیت آماری پارامترهای متغیر 


انحراف 
خی ۰ ۳۳۳ میانگین 
زا بوع نکین فا 
(ه ول | تصادفی نرمال جدول (۳) ۳/۲۸ 
۳ جدول 
(ورثص) ,0۵ | تصادفی | لوگ نرمال | جدول (۲) 5 
۵ تصادفی نرمال ۳/۸ ۱/۰۹۵۹ 


تحلیل فراوانی سیلاب 
برای تخمین حداکثر دبی سیل در دورة بازگشت‌های 
مختلف در تحلیل فراوانی تک متغیره یلاب توزیع‌های 
آماری مختلفی بر مقادیر حداکثر سیلاب سالانه برای یک 
دور 44 ساله (از سال ۱۳4۶ تا ۱۳۹۲) برازش داده 


شده است. توزیع‌های مورد استفاده لوگ نرمال, پیرسون 


٩۱ 


تیپ پنج سه پارامتره لگاریتم پیرسون تیپ 111 لوگ نرمال 
سه پارامتره» پیرسسون تیپ پنج. گامبل حداکثره لوگ - 
لوجستیک و مقداری حدی تعمیم‌یافته هستند. 

نتایج آزمون نکویی برازش با استفاده از نرم‌افزار 
1 ۳8 برای انتخاب توزیع مناسب بر اساس آزمون کای 
مربع در جدول (۵) آمده است. نتیجه آزمون ذشان می‌دهد 
که همه توزیع‌ها برای داده‌های سیل مطلوب هستند و 
هیچ‌گونه اختلاف معناداری بین توزیع تثوری و توزیع 
داده‌های تجربی وجود ندارد اما همان‌طور که مشساهده 
می‌شسود توزیع لوگ نرمال برای سیلاب برازش بهتری 
نسبت به دیگر توزیم‌ها دارد. 

تحلیل فراوانی دو متغیره سیلاب نیز با استفاده از توزیع 
لوگ نرمال صورت گرفت. علت انتخاب توزیع لوگ نرمال 
دی متخیره ابر است که ماس رین« توزیع برای, مر دی 
اوج بوده و بر متغیر حجم سیلاب نیز برازش خوبی یافته 
است. نتایج آزمون نکویی برازش برای پارامتر حجم سیلاب 
نیز بر اساس آزمون کای مربع در جدول (1) ارائه شده است. 


جدول ۵ برازش توزیع‌های مختلف آماری بر حداکثر مقادیر سیلاب سالانه (50*/5) 


توزیع احتمال 


پیرسون تیپ ۷ 
لوگ نرمال 


پیرسون ۷ سه پارامتره 

مقداری حدی تعمیم‌یافته 
لوگ - لوجستیک 

لوگ نرمال سه پارامتره 


لکاریتم پیرسون تیپ 111 


گامبل حداکثر 


آزمون کای مربع 
مقادیر آماری مقادیر نتیجه | رتبه 
و 

۱/۰۳ ۷ | مطلوب | ۷ 
۳-۰۱ ۸۷ مطلوب | ۱ 
2 ۷ | مطلوب | ٩‏ 
1-۹۱۳۹ ۷ | مطلوب | ۲ 
1 ۸۷ مطلوب | ۳ 
۹۹۸ ۷ | مطلوب | ؛ 
۱/۱ ۷ | مطلوب | ۵ 
۱۷/۳۳۳ ۹/۸۷ مطلوب | ۸ 


جدول 1 نتایج آزمون نکویی برازش بر حداکثر مقادیر سیلاب سالانه و حجم سیلاب 


سال سی و دوم شمارة دو ۱۳۹۸ 


حجم سیلاب (۳) 


توزیع لوگ نرمال 
مقادیر آماری مقادیر بحرانی نتیجه 
۱/۱9۰۵ ۱/۷ مطلوب 


نشریه مهندسی عمران فردوسی 


۹۲ 


علاوه بر این شکل (۵) نمودار ۳۳ (- بانلنطادطاظ 
۲ را بر اساس توزیع احتمال اتخاذشده به ترتیب 
برای م00 و ۷ نشان می‌دهد. نمودار ۳-۳ نمودار احتمال- 
ایرال و با فرص فرص اسر حارط نا ادها 
تابع توزیم تجمعی "1 ساخته شده است. اگر روند نمودار 
خطی باشد بدان معنی است که مدل از توزیم در نظرگرفته 


شده پیروی می‌کند. 


۰۰ 
نت 


۰ 

0 

الف) ش ۱ 
0 


01 ۳ 053 0" ۳۷5 13 07 08 059 1 
۱ 


هه عد 


(1001) ۳ 
و هو 
۳ 
ُ 
۰ 
۰ 


1 


1 0/9 05 07 06 0۷5 0۳ 03 02 01 
(ازنوه۲) ۲ 


شکل ه نمايش ۳-۳ بر اساس توزیع لوگ نرمال الف) دبی پیک 


سیلاب ب) حجم سیلاب 


کی ۷ 1-50 ۰ داده های مشاهده ای ۰ 


۷۵ 1-200" - - ۷۵۷۲۵ ۲2100" سب 


36 40 44 48 52 56 60 4 


جدس 1 ۳ کر ۳ 
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ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 


تحلیل فراوانی دو متغیره سیلاب 
تحلیل‌های توأم بر خلاف تحلیل‌های تک متغفیره که تنها یک 
عدد را به ازای هر دورة بازگشت در اختیار قرار می‌دهند. 
بازه‌ای از مقادیر را برای هر متغیر به ازای هر دوره باز گشت 
در اختیار قرار می‌دهند. مجموعه‌ای از دبی اوج و حجم 
سیلاب با دور بازگشت توأم محاسبه شده و در شکل (۷) 
رسم شده است. در این مطالعه از دوره بازگشت توأم "یا" 
برای دبی اوج و حجم سیلاب استفاده شده است. همان‌طور 
که در شکل (1) مشاهده می‌شود. منحنی دورة بازگشت توأم 
می‌تواند به‌طور جانبی در امتداد محورها گسترش یابد 
درساتی که نیز نان آهمیت داهیها: فرمت ای فایل قیرق 
است و باید توسط حدود پایین و بالا محدود شوند. از این‌ری 
بر اساس روش پیشنهادی 16ظ13 [21] منحنی‌ها با حطوط 
عبوری از مبداً با شیب‌های حداکتر (سسمت) 00:۷ و حداقل 
(«نمت) 0۳۷ محدود می‌شوند. همان‌طور که در شکل 
مشخص است ۲۵۶ و من« که به ترتیب شیب حداکثر و 
حداقل نامیده می‌شوند. خطوط ماکزیمم و مینیممی هستند 
که از مان داه‌های مشاهداتی خوو کرناند. رام دفووه 
بازگشت‌های ۸۵۰ ۱۰۰ و ۲۰۰ سال شش نقطه (671-0 تا 
(۷6۵) با مقادیر مشخص حجم و دبی در نظر 


گرفته‌شده‌اند. 


1500 
1350 
120 


1050 


0) 


0 4 8 12 16 20 24 28 


۲ )۲۱۵۲۷( 


شکل " منحنی دور بازگشت توأّم حجم سیلاب- دبی اوج سیلاب 


نشریه مهندسی عمران فردوسی 


سال سی و دوم» شمارة دو ۱۳۹۸ 


آتوسا رضازاده عنبرانی- محمود فغفور مغربی 


حال با استفاده از روش پیش:هادی ۸102772 و 
72 [22] هیدرو گراف‌های مربوط به آن‌ها با استفاده 
از روابط زیر به دست می‌آیند: 

< (۷ ,0,۲ )0۵ 
3 2 
آب ,10 ۱۲ (2)6 - (3)6] مه 
8 8 


روا 2 راعن ۱ 
ف ات تس 02۷01 2(و- 2۷01 ۳ 0 
(0,وت) لا (2-۵۵,0) 6 ]0 


۳۳۱ 
در آن ما و ما به ترئیب زمان اوج و زمان پا یه 


هیدروگراف بوده که می‌توانند به صورت زير محاسبه شوند 


231 

2۷ 
۳ سر - ما 
(۲۵) نز < را 


ُ تحلیا فراوانی دو متغیره دبی اوج و حجم سیپلاب 
به‌دست‌آمده برای دورة باز گشت‌های ۰ ۰ و ۲۰۰ سال 


به ترتیب در شکل (۷) و همچنین جدول (۷) ارائه‌شده‌اند. 


الف) 0 سای 
۷2-- 
دج .- 
4ب 
0 وس 
0 


متوسط هیدروگراف ورودی از تحلیل تک متغیره سس 


35 40 45 50 >> 60 


1 


هکت 
04 
05 وس 
تحت 


متوسط هیدروگراف ورودی از تحلیل تک متغیره سس 


نشریه مهندسی عمران فردوسی 


0 ت 50 که 40 کد 30 


سل 


۹۳ 


ممان‌ظری که ی فزان. ان شکل. (۵ مشاهده کزی 
تفاوت‌های اصلی میان هیدروگراف‌های ۷1-۵ تا ۷6-۵ در 
هر دور بازگشت در دبی اوج خود و زمان پایه 
هیدروگرافنه است, .یرای عال. هیدرو گرات: ۷1-0 که 
کمترین دبی اوج با حجم سیل نسبتاً بزرگ است. در نتیجه با 
زمان پایه طولانی‌تری است» درحالی که هیدروگراف ۷6-۵ 
دارای بالاترین دبی اوج با زمان پایه نسبتاً کوتاه و حجم سیل 
کمتری است. در بخش زیر احتمال روگذری برای همه 
هیدر و گراف تولید شده از تحلیل فراوانی سیلاب تک متغیره 
و دو متغیره و بدترین شرایطی که می‌تواند برای سد تبارک- 
آباد رخ دهد ارائه شده است. 


جدول ۷ دبی اوج سیلاب و حجم سیلاب بر اساس تحلیل فراوانی 


دو متغیره 
(هع) 1 ( ۸ () و ۷0۷ ۷ 

۱۰/۹ -- ۷ ۳۶۵/۵ -- ۶ 0۰ 
۱۳/۲۱ ۲ ۶۳۳/۲ -- ۲ ۳۳ 
۱۵/۳۲ -- ۵۳۳/۱ ۳ ۳.۰ 

1000 

900 

5 هه 900 
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سال سی و دوم» شمارة دو ۱۳۹۸ 
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ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 
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شکل ۷ هیدروگراف ورودی بر اساس تحلیل فراوانی تک متغیره و دو متغیره برای دورهُ باز گشت‌های الف) ۰ ب)۱۰۰ و ج) ۰ سال 


جدول ۸ ریسک روگذری برای ارتفاع‌های مختلف تراز اولیه سطح آب و دورة بازگشت‌های مختلف 


۹۴ 
9 ۰2۷ موه 
2---- 
۷3 ..-- 
۷4 
) ۰ #- 
۷ ۰ 
عتوسط هیدرو گراف ورودی از تحلیل نک متفیره سس 
0 كت "۳ 45 40 23 
دورةٌ روش شبیه- 
بازگشت سازی 2/۱ 13 
٩2-۳۰ ۱/۱۵‏ 0۱۱۰ 
0۰ 
۱+۵ ۰ ۸ ۱22 
ده ۰( ۱۰ 
۱۰۰ 
۱۵ ۰ ۱ ۰ 1-۹( 
ده ۱ ۰ 
۲۰۰ ۰ 
۱*۵ ۱-۰( ۱ 


ریسک روگذری بر اساس تحلیل فراوانی یک- 
متغیره سیلاب 
احتمال روگذری در اثر سیلاب با دورة بازگشت ۰ ۱۰۰ 
و ۲۰۰ سال با در نظر گرفتن سه متغیر غیرقطعی شامل دبی 
اوج, سطح اولیه آب و ضریب تخلیه سرریز محاسبه شده 
است. تمام متغیرهای غیرقطعی, مستقل در نظر گرفته شده 
و برای شبیه‌سازی» روش‌های مونت‌کارلو با حجم نمونه 
۰ و نمونه‌گیری لاتین هایپرکیوب با حجم نمونه 
۰ استفاده شده‌اند. ریسک روگذری در دوره بازگشت 
های مختلف و مقادیر مختلف سطوح اولیه آب برای هر دو 
روش شبیه‌سازی مونت‌کارلو و لاتین هایپرکیوب در جدول 


(۸ اراتد گم سکاو انم ال هی فان یه کرقیق که 


نشریه مهندسی عمران فردوسی 


تراز اولیه سطح آب (0) و 


11 4 ۹ ۸ 
(۱۱ ۰ (۱ ۰ ۱/۱ ۹۰ 0۱ ۰ 
۱۰ (۲ ۱ ۰ ۱ 2 
۱/۳ ۱2۱ ۱۰ ۱ ۰ 
(۱۰ ۱ ۱۰ ۰ 
۱۳۰ ۱2 ۱ (1-۰ 
۱/۹ 32 ۱۵۵ (۱ ۰ 


با بالا رفتن تراز اولیه سطح آب مقدار ریسک شکست سد 
افزایش می‌یابد به‌طوری که در حداکثر ارتفاع 11 متر احتمال 
روگذری بسیار بالاست. همچنین با افزایش دور بازگشت 
مقتان رسک افزارشی بافته سکم 


ریسک روگذری بر اساس تحلیل فراوانی دو 
متغیره سیلاب 
تحلیل فراوانی سیلاب تک متغیره اغلب بر مقدار اوج سیل 
تمرکز دارد و یک ارزیابی محدود از سیل فراهم می‌کند. 
به‌عبارت‌دیگر این روش یک هیدروگراف واحد ایجاد می 
کند که تنها یکی از بسیار هیدروگراف‌های ممکن حوضه سد 


را نشان می‌دهد. از این رو تحلیل دو متغیره سیلاب برای 


سال سی و دوم شمارة دو. ۱۳۹۸ 


آتوسا رضازاده عنبرالی- محمود فغفور مغربی 


ارزیابی توزیع توأّم دبی اوج و حجم سیلاب و در نظر گرفتن 
هیدر و گراف جریان بیشتر برای تجزیه‌وتحلیل ریسک و عدم 
قطعیت به کار گرفته شده است. ریسک روگذری به علت 
سیل‌های مختلف برای دور باز گشت‌های ۸۵۰ ۱۰۰ و ۲۰۰ 
سال در شش سطح اولیه آب. با دو روش عدم قطعیت 
مونت‌کارلو و لاتین هایپرکیوب بررسی شده و نتایج در 
ول اه که شا اد رهگ مامتها و را 
هیدرورگراف‌های 0 - ۷1 تا 0 - 76 دارای روند تغیبرات 
یکسان نیست. به‌عنوان‌مثال در دوره بازگشت ٩۰‏ سال در 
روش مونت‌کارلو ریسک محاسبه شده برای هیدروگراف- 
های 0 - ۱:1 ۵ - ۷5 در همه ارتفاع‌ها به‌جز ارتفاع ۳" 
متر روند نزولی داشته اما در روش لاتين هایپرکیوب ریسک 
تا دای کات ها 0 تس ور 
وروند نزولی داشته است. در نتیجه نمی‌توان قانون کلی 
برای روند تغییرات ریسک در هر روش در نظر گرفت؛ و 
به همین ترتیب در دور بازگشت‌های دیگر نیز روند 
تی رشن کف یریگ آف‌ها معاوت اشران 
نکته قابل ذکر است به دلیل اينکه هیدروگراف 0۱۷- 
دارای بیشترین حجم سیلاب بوده در نتیجه در همه ارتفاع‌ها 
از ریسک بیشتری نسبت به بقیه هیدروگراف‌ها برخوردار 
است. همچنین همان‌طور که در ریسک روگذری بر اساس 
تحلیل فراوانی تک متغیره نیز مشاهده شد در همه دور 
بازگشت‌ها با افزایش ارتفاع مقدار ریسک نیز افزایش یافته 
است و در ارتفاع 11 متر مقدار ریسک روگذری برای سد 
فا 

شکل (۸) روند تغییرات ری سک روگذری در مقابل 
عمق اولیه آب برای دور بازگشست‌های ۸۵۰ ۱۰۰ و۲۰۰ 
سال به ترتیب برای روش مونت‌کارلو برای هیدروگراف 
جر یان تک متغیره(ی0) و دومتغیره 0 - 0.۷1 - ۷2 
و ها کی زا هو کرایت جر بان 
تک متغیره (,0) و هیدروگراف جریان دو متغیره 0 - ۷1 
۷2 رشان می‌دهتن: خمان‌طور ده 


می‌شود با افزایش ارتفاع آب در مخزن مقدار ریسک 


سال سی و دوم شماره دو ۱۳۹۸ 


۹۵ 


افزایش می‌یابد. همچنین هیدرو گراف‌های با حجم رواناب 
بیشتر ( - 71و00 - ۷2) ریسک بیشتری نسبت به دیگر 
هیدروگراف‌های جریان دارند و حطر سرریز شدن آب از 
بدنه سد خاکی در این حالت‌ها بیشتر است؛ بنابراین در هر 
و ادن اه تیوک تفای 
تحلیل فراوانی سیل دو متغیره بیشتر از نتایج تک متفیره 
است. همچنین با افزایش دور بازگشت نیز مقدار ریسک 
افزایش می‌یابد. 

شکل )٩(‏ ریسک روگذری هیدروگراف ورودی-۷1 
۵ را برای دور بازگشت‌های ۵۰ ۱۰۰ و ۲۰۰ سال در 
مقابل ارتفاع اولیه سسطح آب به روش مونت‌کارلو نشان 
می‌دهد. همان‌طور مشاهده می‌شود با افزايش ارتفاع اولیه 
سطح آب و افزايش حجم سیلاب و دبی ورودی به مخزن, 
تک کت مت بر ی یتیک ای ۳0 افانشی 
می‌یابد. 

مک هام( زا هید گرافاه تیان جک 
متغیره و دو متغیره و ریسک روگذری را به ترتیب برای 
هی راهان ار موی ها یی 
بازگشت ۱۰۰ سال نشان می‌دهد. نمودارهای میله‌ای نشان 
دهنده مقایسه ریسک روگذری محاسبه شله هیدروگراف 
سیل ورودی تک متغیره (00) و دو متغیره است. همان‌طور 
که مشاهده می‌شود مقدار ریسک روگذری در تمام ترازهای 
اولیه آب با استفاده از تحلیل فراوانی تک متغیره در هر دو 
روش مونت‌کارلو و لاتين هایپرکیوب کمتر از نتایج حاصل 
از تحلیل فراوانی دو متغیره است که اهمیت آنالیز دو متغیره 
سیلاب برای روگذری سد را نشان می‌دهد. 


نشریه مهندسی عمران فردوسی 


۹۶ ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 


جدول ٩‏ ریسک روگذری بر اساس تحلیل فراوانی دو متغیره برای دور بازگشت‌های مختلف 
دورة روش تراز اولیه سطح آب (0) 1 
بازگشت شبیه‌سازی ۰/۱ 11 5۸ ۱ "1 1 
۶ ۰ 1/۰۵۱ ۱/۳۵۳۱ از ۱۱۵۲۹۱ زا ۳۸۸۰۱۷۱۱ 
1 111 ۲ ۵ | ۱۲۲۱۰۱ از ۱۵۱۹ ۳/۵۱۸۲ 
115 ۱3 ۲/۳۲۹۱ | ۱ | ۵۵۳۱ از ۲۸۵۹۱۰۱ 
ِ 11 ۲۱۲۸۵۱ | ۱۵۵۱ ۵۱۱۱ ۲۸۵۱۷۱۳۰۱۱ 
115 ۱2 ۰( ( ۸۷۷۷۵ ۵۲۱۸۹۱ | ۲/۲۵۲۹ 
۳-50 ٍِ 1115 ۰ ۱/۵۵۱ ۳۸۱۱ زا ۲۱ ۲۸۱۵۹۹۳۱ 
25 ۷/5 ۰ ۱ ۱ 
ِ 1115 ۲/۸۱ ۱۸۱۳۱۳۹۱ ۳۵۳۹ | ۵۲۱۲۵۱ از ۲۸۲۱۷۱۳۴۰۱ 
65 ۱ ۰ ۱/۵۳۵۱ | ۳۲۳۳۱۱ | ۳۸۱ ۱ ۱/۹۹۵۱ 
ِِ 115 ۱2 ها ۱/۵۸۰۱ ۳/۱۱۲۱ | ۳۸۵۱ از ۱/۹۰۵۱ 
۱/۵ ۱۹۰ را ۱/۵۵۲۱ از ۳۲۹۹۵۱ زا ۵۱ ۲۸۰۱۳۰۸۱ 
1 115 ۱ ۰ ۵۵۳۰ | ۳/۲۵۲۰ از ۵۵۵۱ ۲۸۰۲۳۷۹۱ 
۶ ۰ اک اااا ‏ ام ۰( ا ا( ٩-5‏ 
پِ# 1115 ۱۰ ۳۹۰ ۳ 0 ۰ ۱/۸۹۸۱ | 1۱۳۲۵۸۲۱۰ 
1/5 ا ۸ | ۳/۳۹۰۱ | ۱۱۱۳۱۰۱ ۸۱۲۲۵۵۱ | ۳/۲۳۵۱ 
> 11۳5 ۳۸۸۸۵۱ از ۳/۲۱۱۱ ۱۱۳۱۵۱ ۱۱/۵۳۱ | ۳۸۵۰۱۱۳۰۱ 
۱ ۱۰ ۱ ۰ ۱۱ ۲/۶۷۷/۹۵۹ 
1-100 س 11۳5 ۲/۲۸۲۱ | ۲۲۱ ۵۳ | ۵۱۳۲۹ ۲۸۵۷۱۵۱ 
نت4 5 ۰ ۲۰۱۲۸۲۱ | ۳/۵۵۲۱ | ۱۸۱۰۱ ۲۸۰۱۵۳۹۱ 
‌ 11۳5 ۰ ۳۵۵۱ | ۱۸۱۳۵۰۱ از ۳۲۹۲۹ | ۱ | ۲۸۱۱۷۱۳۵۹۰ 
115 ۱2 
۳ 11۳5 ۱ اا ‏ اش( ۱ 
1/5 ۰ ۰ ۲۱۲۸۱۰۱ ۳/۸۵۳۱ | ۵۱۳۱ | ۲/۳۲۵۱ 
ٍِِ« 11 ۰ ۲ ۱۵۱۲۵۰۱ | ۳/۵۵۵۵ | ۵۲۱۹۰۱ ۲/۲۱۸۵ 
۱ ۱/۸۰۰ ۱۳۰ ۲/۸۳۱ | ۲/۲۱۵۰ از 1/۱۵۱۱ 
ِ_ 1115 ا ۷/۱۱۳۱ | ۱۹۱ | ۱۱۲۱۳۲۱ ۱۰۱۵۱ از ۵/۵۱۲۱ 
5 ۰ ۹ ۱ ۰ ۵--* 
تِ 1115 ۵۱۵۱ | ۸۵۱ | ۵۵ ۹ از ۶۷۸۳۷۹۱ 
5 ۱۰ ۰ | ۷/۱۸۵۱ | ۳/۵۵۵۱ 
200 -۳- ِ 115 ۱ 
نت ۶ ۱ 
ِِ« 115 ۸ ۳/۱۵۸۱ ص۵۳۱ ۵ ۳۸۵۱۳۷۰۹۱۱ 
۶ ۲/۷۱ | ۲/۲۱۵۱ | ۵۱۳ ۱ ۲۸۸۰۵۲۰۱ 
۳ 115 اک ا( ۹ اک اا ۲( 
۶ ۸ ۲/۵۱ از ۲۱ ۱۱۲۱۱۱ ۲۸۹۵۱۵۳۱ 
ِ 11 ۱۵۵ 


نشريه مهندسی عمران فردوسی سال سی و دوم. شماره دو» ۱۳۹۸ 


آتوسا رضازاده عنبرانی- محمود فغفور مغربی 
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شکل ۸ تغییرات ریسک روگذری در برابر ارتفاع اولیه سطح آب برای دور بازگشت‌های مختلف 


1 3 5 67 
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25-0 
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شکل ٩‏ ریسک شکست سد هیدروگراف جریان ۷1-0 برای دورةٌ با ز گشت‌های مختلف به روش مونت‌کارلو 


سال سی و دوم. شمارة دی ۱۳۹۸ 


۹۷ 


نشریه مهندسی عمران فردوسی 


۹۸ 


ارزیابی ریسک روگذری سد تبار کآباد بر پایه تحلیل فراوانی یک.. 


56 64 هم و 54 501 
4 ۱ ۲ 1 ۲ ۲ 600 
1 ریسک روگذدری اه تک ریسکت روگذری 2۷/1-0 22 
هیدر و گراف سیل تک متغیره جج یرو کات ۱۷1-5 سح 
10 4.92922 500 
۸.0672 1.19422 
102 ‌_ 


1.0016-3 


۱ 
۱ 1.41664 


وز یآ یر جی: | 


۳ 


40 


(ریل) 2 30 


5.930]-05 


6509-06 


6.153]-7 


عمعجووء ۱ 


10 


شکل ۱۰ ریسک روگذری :0 و ۷1-۵ بر اساس روش ۷105 برای دورهٌ بازگشت ۱۰۰ سال 


56 64 00 ملگ وء هد 501 
10 ۰ ۰ ۰ ۰ + 1200 
1 ریسکت رو گدری ۵۷ ات ریسکت روگذری 226-0 2 
19۰ هید رو گر اف سیل تکت متغیر و << هید رو گر اف 6 - ۱6 سس [ 1000 
2-ع2.325 
عم 7:3:067-02 03-]5.3544 
2.15863 سح 
3 103 3850-04 تفت ‌ 
13 ممععنمد_ حححج 3 
1 10 ومعع2.12 اضع ۷۹ 
۶۰ 10-5 6.50966 سح 1 ‌ِ 
۳۹ ۱ ۵ 
10 7عجمو.د 100 
197 ۱ عمعوصده 1 ل 
۲ 
105 200 
10 تدم 
1092۹ جح < 2-2-2 7.-.۰-.-بسبپبسس__س<. هه 
30 25 20 15 10 5 0 
[ 0 
7 ۳ نز 
شکل ۱۱ ریسک روگذری :0 و ۷6-۵ بر اساس روش ۷05 برای دور باز گشت ۱۰۰ سال 
56 64 (0 ولا مد هد 501 
10 500 
1 ریسکت روگذری له تک ریسکت ر وگذری ۷2-۸ 22 1 وج 
۸ 
10 6-02 3.507 هیدر و گر اف سیل تکك متغیره ----< 1 77 سح ۳۳9 
و رح تحت / 
تست 1 و-ع7.542 1 ۳9 
2589609 ۱ ۲ 1 ۱ 
۶ ۱03 ۳ [ مود چ- 
۳ ۳ 
4 
3 0 05-ع3.426 1 250 ۷4 
3 105 2986و 7 3 
1 200 
" ۱ به 
10 ۱ 150 
1 
107 ۲ ۳ 
1 
1 
10۶ ۱ 50 
10 ت۰۳ ۳۲ ۳ لاسام 


30 


10 


شکل ۱۲ ریسک روگذری :0 و ۷2-0 بر اساس روش 1.15 برای دورة بازگشت ۱۰۰ سال 


بحث و نتیجه گیری 
این مطالعه ریسک روگذری سد خاکی تبارک‌آباد واقع در 


خرا سان ر ضوی را بر | ساس تحلیل فراوانی سیلاب تک 


نشریه مهندسی عمران فردوسی 


متغیره و دو متغیره با تو جه به متغیر های عدم قطءعیت 
مختلف که شامل دبی پیک سیلاب تراز اولیه آب در 


مخزن و ضریب تخلیه سرریز است. ارزیابی کرد. نتایج 


سال سی و دوم» شمارة دو ۱۳۹۸ 


آتوسا رضازاده عنبرالی- محمود فغفور مغربی 


ریسک روگذری بر اساس تحلیل فراوانی تک متغیره نشان 
می‌دهد که هر دو روش مونت‌کارلو و لاتين هایپرکیوب. از 
یک روند یکسان پیروی می‌کنند و ریسک شکست سد با 
بالا رفتن سطح اولیه آب و دبی های ورودی به مخزن 
افزایش می‌یابد. همچنین مقدار ریسک در ترازهای بالاتر 
آنشتی انیتخ: 

تحلیل فراوانی دو متغیره سیلاب با | ستفاده از توزیع 
لوگ نرمال دو متغیره با در نظر گرفتن شش هیدروگراف 
ورودی ۷1-۵ تا ۷6-۵ برای دور بازگشت‌های توأم ۵۰ 
۰ و ۲۰۰ سال انجام شد. هیدروگراف‌های با حجم 
رواناب بیشتر (-۷1 و ۷2-۵ ) ریسک بیشتری نسبت به 
دیگر هیدروگراف‌های جریان دارند و حجم سیلاب بیشتر 
باعث افزایش خطر روگذری می‌شود. نتایجم حاصل از 
روش ۷۵0۳۱۵7۶ و 00 [9] بر اساس تحلیل فراوانی 


۴ 


۹۹ 


دومتغیره یز نشان داد حجم سیلاب کمتر باعث 
اطمینان‌پذیری روگذری سد است. 

همچنین مقادیر محاسبه شده ریسک در تمام ترازهای 
اولیه آب در هر دو روش مونت کارلو و لاتين هاییرکیوب 
با استفاده از تجزیه‌وتحلیل فراوانی سیل دو متغیره بیشتر از 
نتایج تک متغیره است. همان‌طور که مشاهده شد در برآورد 
ریسک روگذری بر اساس تحلیل فراوانی تک متغیره و دو 
متغیره با بالا رفتن تراز آب در مخزن, به‌مراتب از افزایش 
بیشتری نسبت به زمانی که سیلاب‌های با دوره بازگشت‌های 
طولانی‌تر اتفاق می‌افتد. روبرو خواهد بود. 

نتایج تحقیق انجام‌شده توسسط گودرزی [11] که 
ریسک روگذری را توسط تحلیل فراوانی دومتغیره با توزیع 
احتمالاتی گامبل لاجستیک برای دوره باز گشت ۱۰۰ سال 


محاسبه کرد نتایج کار حاضر را مورد تائید قرار می‌دهد. 


۰ 296) م۵۳00 0۷ صحل ۲0۲ ولو امه واصته) همجن مصه لفنک و۷ رتفمل2 20 رآ .سا تاه بط 000022 .1 


2014(۰) ,30-61 .00 و6 ۷۵۰ ,65۵0۲6 ۵۷۱۲۵0۵۱۵- ۵ 0۴ ۵9۲۱۱۵ ,ص1۳2 بجند1 120۲۵۱0۵02 م1" 


۲ سالاری, م. آخوندعلی ع. م. ادیب, آ. و دانشخواه ع. "حلیل فراوانی سیلاب دو متغیره با استفاده از توابع مفصل " علوم و مهندسی 


آبیاری» ۳۷ صص ۰۳۸-۲۹ (۱۳۹۳). 


۰ 2002 ۵۴ 0۵ ولو امه تمصمنوع 11000 لمممنعه: ۵۶ فاتتمجط صرح ولمطاه]۷ بت مقصقصصت 
.(1988) ,269-290 00۰ ,100 

1(۰) ۳۱۷۲۵0۱۵2۷ طا و1۳۵ رومزنهو 11000 تمصصه ۵۶ فلورلممه لمعتاعتاهاگ ر.ظ رصعمعووتصوه۴ م2 .ظ ,ععاظ 
1994(۰) ,117-135 00۰ ,1۱012 رجمت2 ما1 متهعومک متتامملمو ۵۶ اتمصتامن) 

۵ 0۵۱۳۱۵ ,۲۱۵۷۷۹ 11000 ۵۶ عصناممصظ مامه ات۱ رو رعتصفط مضه .و رطاعگ رل بلع00 
1998(۰) ,146-155 .00 ,124 ,۷۵۱ ,۲۱8۵1۱۱66718 

۰ ,1۱۲۵۱۵110021 ماج ۷۷ م فلو مه برمجمنامع 11000 ما ممتاتا ماوت تقصیمه مهتنه 0۱ ما راممرم ریق رمتا لا 
(1999) ,248-254 .00 ,24 

لهمزع1۳010 ۲۱۲ ,-18006و6 11000 مفهزنه 1ص چ عم مه ممتاتاتمافتل لقممصوما ملمتیهژه فط1 رو رما 
2000(۰) ,2575-25868 00۰ ,14 ۷۵۰ ,۲۲0۵6۵8565 

طمومتممره آهعتاعتاهاه ماقتنه17ظ رب رمووم 0ص بظ رقلهع2عظ .۷۲ ,تمه ,.ت ,21720011 , ,۱:6۲6[6 ع(] 
2005(۰) ,50-97 00۰ ,10 ۵1۰ ۷ ۱6۱8و و0۳۵0 0۵ 5۵111۷2 صصعل ۵۴ بدمهنا0ع20 >01عوه ما 
0 2 0۲ وارلمصهه مه متام ممتممماو۲ن ۵۶ اصمصصووموو۸. ر.ظ ۷ رتصتصت مضه ۷۲ م۸ ر2فصصه ۲ 


2008(۰) ,1177-1182 .00 ,35 ۷۵۱۰ ,۲۱۵۱۱6۵8 ۷ ۵۴ 0۱۵ 6۵۳۵۵۵ حصصجل حمتاجعاع 


سال سی و دوم شماره دو ۱۳۹۸ 


نشریه مهندسی عمران فردوسی 


۱۰ ارزیابی ریسک روگذری سد تبارکآباد بر پایه تحلیل فراوانی یک... 


۰ منصوری, ن.» کبیری سامانی» ع. اتحلیل ریسک روگذری سیل در سد با در نظرگرفتن عدم قطعیت‌های محتاف به کمک روش 


شبیه‌سازی مونتکارل و " ششمین کنگره ملی مهندسی عمران. دانشگاه سمنانه (۱۳۹۰. 
0ص ۷۵26 رم 0ععهها کلف عصام‌ممه۲ن حصعل ۵۴ ممتاه2 ۳۲ ال رتفمل2 مصه ,۱۷۲ بتعه۷۱۲2 ر.ظ ق06۵۵ 
,374-387 .00 39 ۷۵ ۱6۵۵و ۱۲ ۵ م0 0۵۵0 و 69و لمصه م60 11000 ماهنته/زه 
2012(۰) 


11۰ 


۲. شرافتی ا. ذهبیون. ب. "تحلیل ری سک روگذری سد با درنظرگرفتن عدم قطعیت‌های هیدرولوژیکی و هیدرولیک ی" مجله علمی - 


پژوهشی هیدرولیک» صص ۸-۱۷ (۱۳۹۳). 

0 09۱۱۵ ,00طامص انامه فمطا عصتمها فتورزلقصه بمممنوع 11000 ماوتنه۱۷ظ. ,۷۰ مهوت مه ما رعصقط2 
2006(۰) ,150-164 00۰ ,11 ۷۵۱۰ ۱61۵و و۱۵0۵ 
رف ۱ ۱ 7 6۱ ( 
(2001) ,1033-1045 00۰ و15 ,۲۷۵۱ ,۲۲۵6۵۵5565 21ع1ع17010 ۲۱۲ 
حصهل طایهه تم فتوولمصه عفن ریق بتقلنتا0 20 ۱۷۲ رعصمعصم رل ملع ما بتاع‌مصه۲ با عمجم 
2008(۰) ,76-87 00۰ ,1 ۷۵۱۰ رعصلمعصنعصظ 0ج ممممزمی تمه ۷۷ فطام۵۵ 0۷۵ 
0ص اصمصوععوعه هناتطمتاه عمتم‌ممتومه فصماورو0 ۳۱۷ رو ی روطتطهام]۷ 20 بب) .ظ روا ویک .۷ ,ع188 1 
2006(۰) ,"20217516 >1و11 
فصح 11 ۶۵۲ اممطصقط 2 :فتورلممه وتلتطمناه: مه عفن ه بتعه1 224 لو رصتهل بیظ.۷ مطعمله 
2007(۰) ,.۵1866۲6ظ 01 0۲ وامه50 مهمتمصد رد وهممتومه لهاممصممتز وم 
۴ ۳۲۵19۱ مان وماز۷۵ن حصقل ۶۵0۲ فلوتاممه عافن۳. ۲۱۰۳۰ بش 0ص2 .۷ رتاقا رن).ظ باه ۷ و.].1 ,ماک 
2007(۰) ,9559-963 00۰ ,133 ,۷۵۱ ,و۱۵۳ و عه رازه ال ,۰ تناو وقهه 2 وه 

.)۱۳۹۲( بزرگ حداد. | مقدمه‌ای بر تحلیل عدم قطعیت در سامانه‌های متاب عآب. تهران انتشارات دانشگاه تهران‎ ٩ 
0600۵022 ما وردم بط ما اصتدانهمصه لمع عافتر مه صمتاملماص1 .۲۲ تتاطو. مصه بآ بلفعل2 و‎ 
ی مممممگ مممتمد رد مصتههمتوجه‎ ۱1510655 1۷۲6012, )2013( 
عطا ۲ ممتتقصصتای مطا ۲۵۲ امعمصمن مفنومل متاونازمهطامنس لمصمزمممصتلل]۳۷ ریم رماطام۲1‎ ۵۲۵۵۵ 
ههانه ۷ ,"عصول ۵۶ انانطهط00م‎ 0. 166۰ ۱۲1۷۰ 0۳22, )2001(۰ 
فص ۶ فعصن0عع۳۵6ظ ,متام ۶۱۵00 صوتععل حصعل ۶۵0۲ طمهمتممج 0۵۷۷ ۸ ره و2عتتصط2؟ 20 بظ رقحصعل‌آخ‎ 
1۸۳۲۳۴ 000821655, بةتاوتاظ و0122‎ 00. 22-27, )1999(۰ 
0002ه) قعمتاموم- مه ۵۶ ات‌ممم‌ممی خر و010 ممناممج گم ا00صاصها ۳ و۷۰۲ ,سمل‎ 
[۷]۵027۷7-۳1111, )1964(۰ 
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